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Abstract 

It is shown that (except for two well defined cases), the necessary 
and sufficient condition for any spherically symmetric distribution of 
fluid to leave the state of equilibrium (or quasi-equilibrium) , is that 
the Weyl tensor changes with respect to its value in the state of equi- 
librium (or quasi-equilibrium). 
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1 Introduction 



Since the publication of Penrose'work |jj , there has been an increasing inter- 
est in the possible role of Weyl tensor (or some function of it) in the evolution 
of self-gravitating systems (0 and references therein). This interest is rein- 
forced by the fact that (at least) for spherically symmetric distributions of 
fluid, the Weyl tensor may be expressed exclusively in terms of the density 
contrast and the local anisotropy of the pressure , which in turn are known 
to affect the fate of gravitational collapse [|J. 

In this work we shall consider spherically symmetric fluid configurations 
which are initially in equilibrium (or quasi-equilibrium), and we shall look 
for necessary and sufficient conditions to leave such regime. Usually, these 
conditions are expressed through the adiabatic index, (or some function of 
it) which indicates the variation of pressure with density, for a given fluid 
element ( see || and references therein). 

As we shall see here, such departure from equilibrium (or quasi-equilibrium) 
is allowed if and only if, the Weyl tensor within the fluid distribution , changes 
with respect to its value in the initial (equilibrium or quasi-equilibrium) 
state. 

There is however two possible exceptions for this. One is represented by 
the "inflationary" equation of state p + p = 0. The other situation when our 
result does not apply, corresponds to the case when then system leaving the 
equilibrium enters into dissipative regime and kT/t(p + p) = 1, where k,T 
and t denote respectively the thermal conduction, the temperature and the 
thermal relaxation time. Both situations imply from the "dynamical" point 
of view, that the effective inertial mass of any fluid element vanish (see || 
and references therein). 

The above mentioned result is reminiscent , in some sense, of a one 
recently obtained , regarding the conditions for a transition from a non- 
dissipative to a dissipative regime in a FRW-flat model O. Indeed in that 
case, the Weyl tensor is always zero, and such transition is only allowed if 
kT/t(p + p) = 1. 

The manuscript is organized as follows. In the next section we give the 
expressions for the field and transport equations, and for the Weyl tensor. In 
Section III the conditions for the departure from equilibrium (and quasiequi- 
librium) are found. Finally in the last Section the results are discussed. 
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2 Field and Transport equations 
2.1 Field equations 

We consider spherically symmetric distributions of collapsing fluid, which for 
sake of completeness we assume to be anisotropic, undergoing dissipation in 
the form of heat flow, bounded by a spherical surface S. 
The line element is given in Schwarzschild-like coordinates by 

ds 2 = e u dt 2 - e x dr 2 - r 2 (d6 2 + sin 2 6d<f) 2 ) (1) 

where v{t,r) and X(t,r) are functions of their arguments. We number the 

coordinates: x° = t; x 1 = r; x 2 = 6; x 3 = 0. 

The metric (jl]) has to satisfy Einstein field equations 

G> = -SttT; (2) 

which in our case read ||: 

_ 8ir » = -I + e-(^--) (3) 

- 8 ^. 1 = 4 + ^(3 + -| (4) 
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-8ttT 01 = -- (6) 
r 

where dots and primes stand for partial differentiation with respect to t and 
r respectively. 

In order to give physical significance to the Tjf components we apply the 
Bondi approach p. 

Thus, following Bondi, let us introduce purely locally Minkowski coordinates 

(r,x,y,z) 
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dr = e u/2 dt dx = e x/2 dr 



dy = rdO dz = rsinOdcf) 



Then, denoting the Minkowski components of the energy tensor by a bar, we 
have 



rpU 



T ° T\ = T\ Tl = Tl T| = T| T 01 = e - {v+x)/2 T { 
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Next, we suppose that when viewed by an observer moving relative to these 
coordinates with velocity u in the radial direction, the physical content of 
space consists of an anisotropic fluid of energy density p, radial pressure P r , 
tangential pressure Pj_ and radial heat flux q. Thus, when viewed by this 
moving observer the covariant tensor in Minkowski coordinates is 



/ p -q \ 

-q P r 

P ± 

V P ± ) 



Then a Lorentz transformation readily shows that 

— p + P r u 2 2Que x / 2 



rpU rpv 

1 o — J o 



T{ = T{ = 
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n = -p± 
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T 01 = e^' 2 T 01 = 



(p + P r )uoe 



{v+\)/2 



Qe _ 

(I-uj 2 ) 1 ' 2 



[l + u 2 ) (10) 



with 



Q = 



qe 



-A/2 



1 -oo 2 ) 1 / 2 



(11) 
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Note that the velocity in the (t, r, 6, <p) system, dr/dt, is related to u by 

w = £ e <*-»0/2 (12) 
dt 

At the outside of the fluid distribution, the spacetime is that of Vaidya, given 
by 

dg2 = ^ _ ^MM\ du 2 + 2dudn _ U 2 ^ 6 2 + sin 2Q d(j) ^ (13) 

where u is a time-like coordinate such that u = constant is (asymptotically) 
a null cone open to the future and 1Z is a null coordinate (#7^ = 0). It 
should be remarked, however, that strictly speaking, the radiation can be 
considered in radial free streaming only at radial infinity. 
The two coordinate systems (t,r,6,<f)) and (u,TZ,6,<p) are related at the 
boundary surface and outside it by 

u = t-r-2Mln(—-l) (14) 

TL = r (15) 

In order to match smoothly the two metrics above on the boundary surface 
f — r T,{t), we have to require the continuity of the first fundamental form 
across that surface. As result of this matching we obtain 

[P r ] s =be^(l-o; 2 ) 1/2 l =[g] s (16) 



s 

expressing the discontinuity of the radial pressure in the presence of heat 
flow, which is a well known result [[|. 

Next, it will be useful to calculate the radial components of the conservation 
law 

= (17) 
After tedious but simple calculations we get 

(-8<)' = ^ (Ti - 7f ) + W (Ti -IJ) + ^(X+|- X -t) (18) 



5 



which in the static case becomes 



K-~(P + Pr) + 2{P \~ Pr) (19) 

representing the generalization of the Tolman-Oppenheimer-Volkof equation 
for anisotropic fluids JjlJ • 



2.2 Transport equations 

As it is well known, the Maxwell-Fourier law for the radiation flux, usually 
assumed in the study of stars interiors, leads to a parabolic equation (diffusion 
equation) which predicts propagation of perturbation with infinite speed (see 



Il]]-||13|| and references therein). This simple fact is at the origin of the 



pathologies |14| found in the approaches of Eckart JT5| and Landau [|l(| for 
relativistic dissipative processes. 

To overcome such difficulties, different relativistic theories with non- vanishing 
relaxation times have been proposed in the past ||T7|-||19|. The important 



point is that all these theories provide a heat transport equation which is 
not of Maxwell- Fourier type but of Cattaneo type pCfl , leading thereby to a 
hyperbolic equation for the propagation of thermal perturbation. 
Accordingly we shall describe the heat transport by means of a relativistic 



Israel-Stewart equation [17] , which reads 



T^ + q a = KP^{Tp-Tap)-ru a q^-\KT 2 (^—u p ) q a (20) 

JJS I \K1 Z / ;/3 

where k, r, T, q 13 and a 13 denote thermal conductivity, thermal relaxation 
time, temperature, the heat flow vector and the components of the four 
acceleration, respectively. Also, P a/3 is the projector onto the hypersurface 
orthogonal to the four velocity u a . 

In our case this equation has only two independent components, which read, 
for a. = 

re {x - u)/2 (Qu + Qu + QuX) +r( Q'u 2 + Quu' + 9^L\ 



r ^ V ) (l-u 2 ) 1/2 (l-cu 2 ) 1/2 
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where the expressions 



g" = Q (^e (A ^ )/2 , 1, 0, 0) (24) 

have been used. 



2.3 The Weyl tensor 

. For the next section we shall need the components of the Weyl tensor. Using 
Maple V, it is found that all non-vanishing components are proportional to 



W= r ~Ci 32 = + + (25) 

where 

^) = ^(p-p + — -T-Y-2^ + 2^j (26) 

corresponds to the contribution in the static (and quasi-static) case . Also, 
the following expression relating the Weyl tensor through the source terms, 
may be found || 

W = -f [ r 3 (T°)'dr + fr* (t| - T, 1 ) (27) 

3 Leaving the equilibrium(quasiequilibrium) 

Let us now consider a spherically symmetric fluid distribution which initially 
may be in either hydrostatic and thermal equilibrium (i.e. uj = Q = 0), or 
slowly evolving and dissipating energy through a radial heat flow vector. 
Before proceeding further with the treatment of our problem, let us clearly 
specify the meaning of "slowly evolving". That means that our sphere 
changes on a time scale which is very large as compared to the typical time 
in which it reacts on a slight perturbation of hydrostatic equilibrium. This 
typical time is called hydrostatic time scale. Thus a slowly evolving system 
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is always in hydrostatic equilibrium (very close to), and its evolution may be 
regarded as a sequence of static models linked by @. 

As we mentioned before, this assumption is very sensible, since the hydro- 
static time scale is usually very small. 

Thus, it is of the order of 27 minutes for the sun, 4.5 seconds for a white 
dwarf and 10 -4 seconds for a neutron star of one solar mass and 10 Km 
radius ||. 

In terms of u and metric functions, slow evolution means that the radial ve- 
locity uj measured by the Minkowski observer, as well as metric time deriva- 
tives are so small that their products and second order time derivatives may 
be neglected (an invariant characterization of slow evolution may be found 
in |2J). 
Thus gl[ 

v \v \ 2 ^ v 2 ^ u 2 ^ u = Q (28) 

As it follows from @ and fllPl), Q is of the order 0(u). Thus in the slowly 
evolving regime, relaxation terms may be neglected and (P0| ) becomes the 
usual Landau-Eckart transport equation . 

Then, using ( p8|) and ([Tj|) we obtain (|19|), which as mentioned before is the 
equation of hydrostatic equilibrium for an anisotropic fluid. This is in agree- 
ment with what was mentioned above, in the sense that a slowly evolving 
system is in hydrostatic equilibrium. 

Let us now return to our problem. Before perturbation, the two possible 
initial states of our system are characterized by: 

1. Static 

u = Q = U j = Q = (29) 

2. Slowly evolving 

Cj = Q = (30) 
Q « 0(u) ^ (small) (31) 

where the meaning of "small" is given by (|28|). 

Let us now assume that our system is submitted to perturbations which force 
it to depart from hydrostatic equilibrium but keeping the spherical symmetry. 
We shall study the perturbed system on a time scale which is small as com- 
pared to the thermal adjustment time. 
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Then, immediately after perturbation ("immediately" understood in the 
sense above), we have for the first initial condition (static) 



u = Q = (32) 

oj « Q ^ {small) (33) 
whereas for the second initial condition (slowly evolving) 

Q « 0(w) ^ («ma«) (34) 

Q^u^O (small) (35) 

As we shall see below, both initial conditions lead to the same final equations. 
Let us now write explicitly eq. (|18D . With the help of (|7|)-(|T0|), we find after 
long but trivial calculations 



PL p'u) 2 2uu)'p 2uu'P r 

' + -7- n + Z + 



l-UO 2 l-UO 2 1-U 2 (l-LU 2 ) 2 

2u 3 m'p 2Q'me x / 2 2Qm'e x / 2 2Quj 2 uj'e x ' 2 
+ (1 - u 2 ) 2 + (1 - u 2 ) 1 ' 2 + (1 - u 2 ) 1 ' 2 + (1 - u 2 f 2 

2 r 47rr 3 , 9 n Que x / 2 12nr 3 ( Que x l 2 \ 2 

r L r-2m v / (1-^2)3/2 r -2m \(\ - oj 2 ) 1 ) 

f , p * u 2 ,, p p v , 2g^eV 2 , (p + gj 1 + ^ 2 m 

+ P + 3 J + Pr ~ P± + 7; ^172 + o 1 2 ^T" 

1 — ur (1 — c<j 2 ) 7 2 1 — or r — 2m 

Quje x/2 m 2vrr 3 / _ 2 x . _ , 1 + u; 2 

+ — 7T7 1 P r + pu; 2 (p + P r .) ? 

(1-cj 2 ) 1/2 r-2m r-2m V (1 - w 2 ) 2 

87rr3 /r» 2\ <5^e A / 2 4vrr 3 „ A/2 , nN 1 + cj 2 

e ~'* + £-^ (36) 



8vrr V 2 2 / 
which, when evaluated immediately after perturbation, reduces to 
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(P + Pr)m + 4vrr , + p2 ^ + 2 (P r - P ± ) = 
r 2 {l — 2m/r) (1 — 2m/r) ^ r T ' r 8%r 



for both initial states. 

On the other hand, an expression for A may be obtained by taking the time 
derivative of (0) 

8 6 ^ p r > i - ^2 2 (l-o; 2 ) 1/2 2 
(p + P r )u X A/2 1 + u 2 x+ (. + p\_" 

+ i-^ 2 2 +Lie a-u,*) 1 '* V r} i 



1 A 1/9 1 

+ (p + P r )u;- z-i + Qe x/2 



uJ 2 



;i-^ 2 ) 2 (l-o; 2 ) 1/2 



+ ^fiz£>, (38 ) 

which, in its turn, when evaluated after perturbation, reads 

A = -87rre^ +A )/ 2 Up + P r ) cj + Qe A/2 l (39) 



replacing A by fl3~9"|) en (|37|), we obtain 

- e (l/ ~ A)/2 P = (p + P r ) + Qe A/2 (40) 
where R denotes the left-hand side of the TOV equation, i.e. 



R 



dP r AnrP 2 
dr 



+ 



P r m 



1 — 2m/r r 2 (l — 2m/r) 
pm 2 (P ± - P r ) 



AmrpP r 
L — 2m /r 



+ 



2m/r) 



r 



PL 



^(P + P r )-^(P±-P,) 



(41) 



The physical meaning of R is clearly inferred from (fTI])- It represents the 
total force (gravitational + pressure gradient + anisotropic term) acting on 
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a given fluid element. Obviously, R > 0/R < means that the total force is 
directed inward /outward of the sphere. 

Let us now turn back to thermal conduction equation (PD|). Evaluating its 
i-component (given by Eq. (pT|) ) immediately after perturbation, we obtain 
for the first initial configuration (static), an identity. Whereas the second 
case (slowly evolving) leads to 

00 (V + Tjj=0 (42) 

which is to be expected, since before perturbation, in the slowly evolving 
regime, we have according to Eckart-Landau (valid in this regime) 

Q = -ne~ x (t> + ^pj (43) 

Therefore, the quantity in bracket is of order Q. Then immediately after 
perturbation this quantity is still of order 0(uj), which implies ((42|). 
The corresponding evolution of the r-component of the equation (^) yields, 
for the initially static configuration 

T Qe x/2 = —kTuj (44) 
where the fact has been used that after perturbation 

Q = =► T' = - T -^- (45) 

For the second case, the r-component of heat transport equation yields also 
([BP, since after perturbation the value of Q is still given by ([|3|), up to 0{uj) 
terms. 

Finally, combining (|40|) and ([14|) we obtain 

e {u ~ X)/2 R 1 

" = — ( 46 ) 



(P + Pr) (l- 



r(p+P r ) 

or 



- e ("- x)/2 R= (p + P r )u(l -a) (47) 

with a defined by 
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Let us first consider the a = case. Then, (JT7|) has the obvious "Newtonian" 
form 

Force = mass x acceleration 
since, as it is well known, (p + P r ) represents the inertial mass density and 
by "acceleration" we mean the time derivative of uo. In this case (a = 0), 
an outward /inward acceleration (Co > 0/Co < 0) is associated with an 
outwardly /inwardly (R < 0/R > 0) directed total force (as one expects!). 
However, in the general case (a ^ 0) the situation becomes quite different. 
Indeed, the inertial mass term is now multiplied by (1 — a), so that if a = 1, 
we obtain that Co ^ even though R = 0. This decreasing of the "effective 
inertial mass density" ( vanishing when a = 1), has been shown to occur in 
the general (non-spherical symmetric case) (see || and references therein) 
Next, evaluating (p5|) immediatly after the system leaves the equilibrium (or 



quasi-equilibrium) and using (|39D and (0), ([47]) m &y be written as 



W - W (s) 

- 3e ( - A)/2 L - 4 UJ ={p+ P r ) Co {I -a) (49) 



27rr 4 

Therefore, Co ^ implies W ^ Wr s ), and viceversa, unless either (p+P,, 
(in the non-dissipative case) or a = 1 (in the dissipative case),. 
Next, assuming uj = 0, it can be shown without difficulty, that 



W-W, 



for any n > 1. Thus, unless either (p + P r ) = (in the non-dissipative 
case) or a = 1 (in the dissipative case), if the Weyl tensor does not change 
with respect to its value in the equilibrium (or quasi-equilibrium ) state, the 
system will not abandon such state. 



4 Conclusions 

We have seen that any departure from equilibrium (or quasi-equilibrium) 
of a spherically symmetric distribution of matter, is tightly controlled by 
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changes in the Weyl tensor with respect to its value in equilibrium (or quasi- 
equilibrium) , as indicated by (f49|) and (loTf ). As an obvious consequence of 
this, it follows that a conformally flat fluid distribution in state of equilibrium 
(or quasi-equilibrium) will depart from such state, only if it ceases to be 
conformally flat (except for the two cases already mentioned). 
Finally, it is worth mentioning that the last term in (^0]) is frequently omitted 
(the so-called "truncated" theory) ||22|| . In the context of this work both 
components of this term vanish and therefore all results found above are 
independent of the adopted theory (Israel-Stewart or truncated). 
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